Fluctuating Mobility Generation and Transport in Glasses 
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In the context of the random first order transition theory we use an extended mode coupling 
theory of the glass transition that includes activated events to account for spatiotemporal structures 
in aging and rejuvenating glasses. We numerically solve fluctuating dynamical equations for mobility 
and flctive temperature flelds which capture both mobility generation through activated events and 
facilitation effects. Upon rejuvenating, a source of high mobility at a glass surface initiates a growth 
front of mobility which propagates into the unstable low mobility region. The speed of the front is 
consistent with experiments on the rejuvenation of stable glasses, which "melt" from their surface. 
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On human length and time scales glasses seem 
static, yet glasses constantly move at molecular scales. 
Molecules in glasses change their locations through rare, 
activated events at a rate which varies throughout the 
glass owing to its aperiodic structure. Motions in one lo- 
cation cause or relieve constraints thereby changing the 
rate of transitions in neighboring regions. From these 
considerations it follows that the apparent mobility of 
regions fluctuates both in space and time and appears to 
be locally transported. In this paper, we explore, using 
numerical methods, fluctuating mobility generation and 
transport in glasses focussing on rejuvenation after heat- 
ing. The framework used is based on the random first 
order transition (RFOT) theory of glasses [11-0] ■ 

The static appearance of glasses has been formally ad- 
dressed by mean field theories such as the explicitly dy- 
namical mode coupling theory (MCT) [^, Q as well as by 
outwardly static approaches [l[ Q • The similarities be- 
tween these different approaches lead to the suggestion 
that they are approximately equivalent Q. These dis- 
parate analogous approaches indeed do provide exactly 
equivalent results for a large class of spin glass models 
with infinite range interactions [1, ■ In mean field 

models, a dynamic transition to immobility occurs dis- 
continuously at a high dynamic transition temperature, 
T(i, below which an exponentially large number of frozen 
states emerges. At a lower temperature, Tr-, the con- 
figurational entropy of these states vanishes and a true 
thermodynamic transition ensues For systems with 
finite range interactions, the transition at Td is wiped out 
by local activated transitions between the frozen states 
RFOT theory argues that activated transitions re- 
sembling nucleation events become the source of mobility 
in deeply supercooled liquids. These thermally activated 
motions would cease at the thermodynamic transition if 
the entropy crisis occurs or otherwise only at absolute 
zero. 

While the activated mechanism of mobility generation 
is a key aspect of RFOT theory, owing to its connection 
with the mode coupling theory, RFOT theory shares with 
kinetically constrained models [lH the idea that mobility 



can be transported i.e. the notion of facilitation. Mo- 
bility transport effect is an essential feature of RFOT 
theory, diminishing the effects of the instantaneous het- 
erogeneity in activation barriers which otherwise would 
lead to a broader distribution of relaxation times than 
observed 

The idea that the mobility of a region of the fluid de- 
pends on the mobility of its neighbors provides the his- 
toric basis for MCT. Bhattacharyya, Bagchi and Wolynes 
(BBW) suggested a way to bridge the gap between MCT 
which is adequate above the dynamical transition and 
mobility generation through activated transitions which 
occurs below the mean field dynamic transition by in- 
cludingan activated event vertex in the self-consistent 



MCT [IJ, UJ. BBW's treatment did not make explicit 
the spatiotemporal nature of the dynamical heterogene- 
ity. To make clearer the connection we can allow the 
mobility to be an explicitly varying dynamical variable 
in the BBW analysis. The resulting coupled inhomoge- 
neous nonlinear self-consistent differential-integral equa- 
tions are cumbersome to study but assuming slow vari- 
ations in time and space in the mobility allows a con- 
tinuum approximation to be made This gives a 
simple description of a rich set of phenomena especially in 
the case of "rejuvenating" glasses that are heated after 
falling out of equilibrium. Wolynes argued on the ba- 
sis of the continuum theory that mobility transport and 
generation conspire to set up moving fronts in heated 
glasses that resemble flames [iBl- Ignition by mobility 
generation occurs at points within the bulk but coopera- 
tive rearrangements occur with especially large rapidity 
at a free interface which acts as an external source for 
a rejuvenation front [l^ . The theory thus predicts then 
there are both heterogeneous and homogeneous mecha- 
nisms by which glass "melts" into an equilibrated fluid. 
These two mechanisms have been observed by Ediger 
in his studies of ultrastable glasses prepared by vapor 
deposition 17[. The approximate continuum treatment 



of rejuvenation agrees qualitatively with the more pre- 
cise later measurements [l7l - [l9j | but according to Ref. [H 
the simple analytical result for the front speed is quan- 
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titatively off by about two orders of magnitude. In the 
first set of measurements 17[ there was also evidence for 
a modest overshoot phenomenon in which, just after a 
front passes, the mobility seems to be higher than it will 
eventually become in the fully equilibrated liquid. This 
overshoot, not anticipated by Wolynes's calculation, is, 
however, not apparent in later experimental studies |18| . 
The possibility of a large deviation from the simple ap- 
proximation inspired us to explore more completely the 
continuum equations for mobility generation and trans- 
port in glasses. Here we show that a more accurate nu- 
merical treatment of the fluctuating continuum equations 
gives better quantitative agreement with observed front 
propagation speed but no sign of the overshoot occurs in 
the average behavior. 

We first recount the continuum equations emphasizing 
the role of fiuctuations. Starting from the spatially inho- 
mogeneous MCT with activated events we approximate 
the memory kernel as being explicitly time and space 
dependent but treat it as locally frequency independent 
thus representing a mobility field ii{x, t). Expanding the 
inhomogeneous equation in Taylor series yields a contin- 
uum equation which we can supplement with local ran- 
dom terms UM 
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where /2 is the uniform solution of the homogeneous 
MCT equation and ^ is the size of a cooperatively 
rearranging region (CRR). We will refer to Eq. ([T]) 
as linearized inhomogeneous mode coupling theory (Li- 
IMCT). In deeply glassy systems jl the mobility from 
the uniform solution of the BBW equations can be rea- 
sonably approximated by the traditional aging formula- 
tion of Narayanaswamy-Moynihan and Tool in term of 
a fictive temperature Tp, which measures the local en- 
ergy, as in RFOT theory by Lubchenko and Wolynes [6] 

ftiTp, T) = ^0 exp |-f^ - ^ fc~ Tp ' }■ The parameter x 
is predicted by RFOT theory to be inversely proportional 
to the heat capacity discontinuity Acp. To complete the 
description we take the local Active temperature to ap- 
proach the ambient temperature through an ultralocal 
relaxation law 
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Since activated motions occurring right at the sur- 
face feel no mismatch energy penalty on the interface 
along the free surface they have to overcome only the 
mismatch energy on the hemisphere facing the bulk 
[l6j . Free surfaces thus represent a source of inho- 
mogeneous mobility. In the bulk itself there are dy- 
namic inhomogeneities from fluctuations in the local fic- 
tive temperature and from the activated events which 
transport the mobility field. We treat the random 



force terms as coarse-grained white noises with correla- 
tion lengths reflecting the length scale of the activated 
events. The intensities of the noises may be found by 
requiring the linearized equations to satisfy locally the 
fluctuation-dissipation relations. The local fictive tem- 
perature fluctuation Si] is taken to be {r]{x,t)r]{x' ,t')) = 
Ac Ivt ~ ~ where TV* is the number 

of molecular units in a CRR. The fluctuations in mobil- 
ity generation 5g and transport 5j arise due to fluctua- 
tions in free energy barrier heights which are the cause 
of stretched exponential relaxation with a bare stretch- 
ing parameter /3o = 1 / a/I -I- {SFi / ksT)'^ 5] . Linearizing 
the equations and treating p, as constant, the fluctuation- 
dissipation relation yields a mobility generating noise 

{Sg{x, t)6g{x\ t')) = 2flfi^ ( - l) ' S{x - x')6it - 
t') and a mobility random flux {Sj{x,t)dj{x' ,t')) — 

- i) (^) ' - x')5{t -t')^. 

While Li-IMCT is a natural perturbation result, it is 
worth noting that to the same order of perturbation the 
rate of relaxation of ji to its local value could be taken as 
jjL itself rather than /I. This leads to a nonlinear transport 
equation f = |, - m (m - m) + ^3 + ^. This 

equation has the structure of the Fisher-Kolmogorov- 
Petrovsky-Piskounov equation [2l], [S^l so we call this 
set of coupled equations the FKPP inhomogeneous mode 
coupling theory (FKPP-IMCT). This equation turns out 
to do a somewhat better job of describing the experi- 
ments than the Li-IMCT does. 

Eq. ([I} resembles the diffusion equation for some kind 
of excitation. Of course, such "excitations" are quite fic- 
tional - there is no regular structure in a glass that can 
be used to a priori define a defect. The formal analogy 
between inhomogeneous mode coupling and excitation 
models explains why such models while being microscop- 
ically unrealistic do manifest phenomena reminiscent of 
molecular glasses. 

Exact analytical treatment of the equations for mobil- 
ity generation and transport is difficult. It has been ar- 
gued, however, that for aging glasses i.e. systems that 
have been quickly quenched to low temperatures, the 
aging dynamics that can be easily studied in the labo- 
ratory will not be very different from what occurs ini- 
tially and so involves primarily the mobility generation 
events themselves [l^. When a glass is heated however 
the autocatalytic interplay between mobility generation 
and transport qualitatively changes the picture. Rec- 
ognizing the analogy to combustion, Wolynes invoked 
some of the simpler approximations used in the analy- 
sis of flames [i^ which predict the front propagates at a 
speed Y^2/3^^^'sh^ ThXs upper bound to the speed de- 
pends only on the mobility in the high T equilibrated 
material, /i*^'^^, but is independent of the mobility of the 
initially prepared low mobility material. Here we test 
that idea by numerical methods. 
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The numerical calculations were carried out using pa- 
rameters referenced to tris(naplitliyl)benzene or TNB 
17| whose glass transition temperature for r = 100 s, 



logZ) [cm^/s] 
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Tg is 347 K, with ^ = 2.5 nm. In ID, the calculations 
use a lattice of 200 cells with a spacing Aa; = 1 nm. For 
the 2D case, a square two-dimensional lattice of 20 x 100 
cells with mesh size Ax — 1 nm was employed. Peri- 
odic boundary condition are assumed for all cases but 
the surface is represented by a region with intrinsically 
high mobilit y g iven by the Stevenson- Wolynes formula 
16|. Owing to the stiffness of the equations 



we cannot study them at as low a fictive temperature as 
the experiment. 

We first report numerical results obtained from ID de- 
terministic theory. Due to the symmetry of the determin- 
istic problem, these also apply in two dimensions. In the 
Li-IMCT case, the interface width as measured by full- 
width-at-half-maximum of the source term — — A) is 
about 5 nm. This scale is only a bit larger than the size 
of the entropic droplets that generate mobility Q . The 
width of the interface from the FKPP-IMCT is as wide 
as 10 nm, more consistent with the experimental mea- 
surement 

The speeds of front propagation at various tempera- 
tures predicted from the ID deterministic Li-IMCT and 
FKPP-IMCT are shown in Fig. [T] Increasing the am- 
bient temperature increases the speed. The predicted 
linear trend of log(wg, ) with temperature T agrees with 
observations. No calculated point differs from the experi- 
mental line [18| by more than a factor of 2. The predicted 
speed from the deterministic Li-IMCT however does vary 
somewhat more with the initial fictive temperature in the 
simulation than in experiment. Wolynes's approximate 
solution on the other hand showed independence of the 
speed on the initial fictive temperature in agreement with 
observation. The discrepancy in the numerical study con- 
trasting with the correct prediction from the approxi- 
mation may arise because the initial fictive temperature 
in the study is not low enough thus leading to a "cold 
boundary difficulty" , a problem extensively discussed in 
combustion theory using intermediate asymptotics [23| . 

We now consider the stochastic models in one and two 
dimensions. While the three dimensional situation is nu- 
merically too demanding for us to treat presently the 
lower dimensional cases already allow a good compari- 
son with experiment. A realization of the front evolu- 
tion from the 2D stochastic version of the Li-IMCT are 
shown in Fig. [2] through the mobility field ji. Initially the 
sample is prepared in an ultrastable phase by simulating 
the stochastic equations for a uniform system at a low 
temperature, T™'*"'' = 310 K. The system is then in- 
stantaneously heated to a higher ambient temperature, 
T = 342 making the former low mobility state un- 
stable. As the fictive temperature increases, the front 
propagates into the bulk leaving the stable state behind. 
Due to the periodic boundary conditions fronts initiate 
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FIG. 1. (Color online) The speed of the mobility front 
from the dynamics of the ID deterministic Li-IMCT (solid 
red markers), the ID deterministic FKPP-IMCT (opened 
red markers), the ID stochastic Li-IMCT (solid green mark- 
ers), the ID stochastic FKPP-IMCT (opened green mark- 
ers), the 2D stochastic Li-IMCT (solid blue markers), and 
the 2D stochastic FKPP-IMCT (opened blue markers). Dif- 
ferent kinds of symbols indicate different initial fictive tem- 
peratures, the triangles are for = ilQK, circles for 
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315K, and squares for T" 



320is:. The ex- 



perimental measurements as reported in Ref. [13 are shown as 
magenta stars and the black stars represent the more recent 
results from Ref. [T^. All lines are drawn as guides for the eye. 
Inset: The front speed versus relaxa tion time. The continuum 
approximation of the front speed \/2/3^^'^'^'^ is drawn as the 
black solid line. Notice the front speed does not scale quite 
as strongly as t^^- 



at both ends and travel towards the "inner" bulk. 

In addition to the fronts emanating from the surfaces 
at high mobility reconfiguration events take place at some 
bulk locations faster than they do at others. Once these 
initially initiated nucleated regions reconfigure, they cat- 
alyze rearrangement of their low mobility neighbors and 
lead to fronts emanating radially from the rejuvenated 
centers - as described earlier [15|. Fast regions cannot 
as easily jump over slow regions which are more stable, 
so they often detour to avoid an immobile region or be- 
come jammed up at the interface between the two re- 
gions. This latter effect in a specific realization can give 
the appearance of the overshoot of the mobility behind 
the propagating front as shown in Fig. [21 Both stochastic 
Li-IMCT and FKPP-IMCT equations give similar fluc- 
tuating overshoot phenomena. 

To assign a typical speed to the stochastic front, the 
front position may be defined using a one dimensional 
integral of a threshold function z{t) = \ dxQ[^{x, t) — 
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] , where B[/(x)] is Heaviside step function, 
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tion barriers [2^ . Shear banding may arise from the 
couphng of stress to inhomogeneously propagating inter- 
nal mobility fronts. 

Financial support by the D.R. Bullard- Welch Chair at 
Rice University and a Royal Thai Government Scholar- 
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FIG. 2. (Color online) Mobility field snapshot at t = 4.6 x lO" 
s. The upper plot shows the mobility field from 2D stochastic 
Li-IMCT. The color scheme denotes the mobility on a log- 
scale with corresponding color bars on the right. The figure 
in the lower panel show average values of the mobility fields 
along the j/-axis. 
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/i""'""""'" is a constant threshold [2J|. The integration 
measures the amount of sample which has not yet been 
transformed much as in calorimetry 19]. Physically z{t) 
gives the depth of the front from each surface. The 
growth front speed is determined hy v = Az/At, where 
At is somewhat larger than a typical relaxation time 
Averaging this quantity over initial times and initial en- 
semble yields the average front speed Vgr- 

The front speeds for ultrastable TNB glasses prepared 
at various initial temperatures from both ID and 2D of 
stochastic Li-IMCT and FKPP-IMCT are displayed in 
Fig. [TJ The resulting speeds do slightly depend on the 
initial temperatures, but depend exponentially on the 
ambient temperature rapidly rising as the ambient tem- 
perature increases. The front speeds for all stochastic 
cases are larger than for the deterministic cases by a fac- 
tor of two. These speeds from the stochastic simulations 
for both dynamics now deviate by less than 50% from the 
front speeds reported in Ref. [iSl. The growth front ve- 
locities reported in Ref. |l7l are about three times greater 
than those found in Ref. [iSj in which it was suggested the 
samples prepared earlier may contain impurities such as 
the vacuum pump oil which increase the mobility. 

In the present paper we have numerically solved fluctu- 
ating continuum equations coming from inhomogeneous 
MCT that were previously introduced to account for spa- 
tiotemporal structures in aging and rejuvenating glasses. 
Numerically accounting for fluctuations in mobility field 
and Active temperature gives closer quantitative agree- 
ment with experimental measurements than the previous 
approximate treatment. 

The front propagation phenomena discussed here are 
of a very general nature. Any external influence that can 
increase mobility, if it is inhomogeneous in space can lead 
to such fronts. For example plasticizer diluents which 
lower Tg will give fronts which will propagate into the 
bulk. The coupled diffusion may give an overshoot phe- 
nomenon like that seen in some of the experiments [25|. 
Stress also acts to increase mobility by lowering activa- 
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